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which include the resummation of soft logarithms and Coulomb singularities through next-to- 
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mt = 173.3 GeV, we find 



-Tevatron 



7 99+0.31+0.71 
■ '-^^-0.47-0.55 



pb 



^LHC . 



162.6+^-4+i5-4 



-7.5- 



13.4 „], 

14.7 



in good agreement with the latest experimental measurements. 
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1. Introduction 

The total top-pair production cross section On has been measured at Tevatron with an accuracy 
Aatt/Ott of about ±7% [|l|, ^ and the two LHC experiments have already reached similar sensitivity 
with the accuracy of some analyses already at the ±6.5% level. With more statistics being 
collected, the experimental error is bound to reduce even further, opening a new era of precision 
top phenomenology. The total cross section, in particular, can be used to measure the top-quark 
pole mass m, in a theoretically clean way, to constrain new physics and to extract information 
on the gluon distribution function (PDF) of the proton. This clearly requires a precise theoretical 
understanding of the ^-production dynamics; more specifically, predictions beyond next-to-leading 
order (NLO) in standard fixed-order perturbative QCD are necessary. 

Near the partonic production threshold, j8 ~ 0, where j3 = 1 — 4mf / s is the velocity of the 
two top quarks, the partonic cross sections dpp'^ttx are enhanced due to suppression of soft-gluon 
emission and exchange of potential (Coulomb) gluons between the non-relativistic top and anti- 
top. These two effects give rise to singular terms, at all orders in perturbation theory, of the form 
a^ln^'^ j8 and a^/jS, respectively. While the hadronic cross section, 

<^tt{s)= £ f ^ dTLpp,{T,^f)dppi^ttx{^s,^f), (1.1) 

where Lppi are parton luminosity functions, receives contributions from regions where j3 is not 
necessarily small, especially for LHC centre-of-mass energies, the region with j3 < 0.3 still gives 
a sizeable contribution to Otj, due to the rapid fall off of the parton luminosity functions at large T. 
One therefore argues that an all-order resummation of soft and Coulomb corrections provides more 
accurate predictions than the fixed-order calculation. 

Leading logarithmic (LL) and next-to-leading logarithmic (NLL) resummation in the so-called 
Mellin-space formalism has been known for a while [^]. More recently, thanks also to the calcu- 
lation of the relevant soft anomalous dimensions ||^, next-to-next-to-leading logarithmic (NNLL) 
resummations [Q, H, ^, and approximated NNLO cross sections constructed from the re-expansion 
of the resummed result [|^, H , 1^, |l3]], have become available. Of the aforementioned works, only 



Ref. [^] provides a simultaneous resummation of soft and Coulomb corrections, obtained through a 



general formalism [ |14| ] based on the factorization of soft and Coulomb effects in the context of soft- 
coUinear effective theory (SCET) and potential non-relativistic QCD (PNRQCD). The formalism, 
and the results of [||], will be reviewed in the following. 

2. Factorization and resummation of the tt total cross section 

The basis for resummation is the factorization of the partonic cross section into short-distance 
contributions, related to physics at the hard scale ~ m,, and effects associated with soft-gluon 



emission, which naturally "live" at a much smaller scale ~ /M,j8 . In [14] it was shown that, at 



threshold, the partonic cross section in fact factorizes into three different contributions, 

app,^,jx{s,^)= £ //;^y(m„At) fd(oJR{E-^)Wf{(o,n). (2.1) 
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The hard function H^p, encodes the model-specific short-distance effects, the Coulomb function Jr 
describes the internal evolution of the tt pair, driven by Coulomb exchange, and the soft function 
W/^ contains soft-gluon contributions. Here R denotes the irreducible colour representation, either 
singlet (R = 1) or octet (R = 8), of the tt pair. 



In the approach adopted here [14, 15 1, the hard and soft functions, Hj^^, and W/^, are resummed 
through renormalization-group (RG) evolution equations directly in momentum space, contrary 
to the conventional formalism, where resummation is performed in Mellin-moment space. The 
relevant RG equations and their solutions were derived, for the general case of massive particle 
pairs HH' in arbitrary colour representations /?,/?', in [14 1. For the soft function W/^ the solution 
to the evolution equation reads 



■ exp [-45(Ai. , M ) + 2a^,; (m. , M )] (^77 , Kv ) ^ ( 



277 



d{(0) 



(2.2) 



In (2^), Sf represents the Laplace transform of the fixed-order soft function W/^. The function 
S controls the resummation of double logarithms, while and 77 resum single logarithms. To 
obtain NNLL accuracy, these functions must be included at the three-loop (S) or two-loop (a^ • and 
77) order, while the fixed-order soft function Sf is required at one loop. An expression analogous to 



(p^) can be derived for the hard function H^p,. 



Resummation of Coulomb effects has been extensively studied in the context of PNRQCD and 
quarkonia physics. The potential function Jr is related to the Green's function of the Schrodinger 
operator -V^/rrit - {-Dr) a,/r [1 + ^(a,)]. 



jRiE) = 2Im G[^|jj(0,0;£) A„c(£) + G^';(0,0;£) + . . . 



(1) 



(2.3) 



In Eq. ( |2.3[ ), G^^^j denotes the LO Coulomb Green's function. 



G. 



(0) 
Cfi 



{-Dr)0Cs 
2^-E/m, 



(2.4) 



Gn'o is the correction to the Green function from the a? terms in 



with £■ = \/f — 2m, ~ m,j8^. G|i]j 
the Coulomb potential, and Anc = I + iff{ccj\nP) represents non-Coulomb contributions which 
enter the cross section first at NNLL order [|l0|]. Notice that for E <0, i.e. below the production 
threshold, Eq. ( pj| ) contains a series of ti bound-state resonances when Dr < 0, which are included 
in the NNLL results presented in Section ^. 

It is often convenient to re-expand the resummed results to construct higher-order approxima- 
tions at fixed order in a^. In particular, at NNLL all singular terms in the hmit j3 — ^ at 0{aj) can 
be correctly predicted. Hence, one can define an approximated NNLO prediction as 



. NNLOap 

V 



y/^NLO , (0) y ^{2,n).nAvl 



(2.5) 



with O'^^ the exact colour-separated NLO cross sections [|l^], and o^^ ^ the Born contributions. 
The NNLO functions /ppf^ incorporate exactly all terms (and only those) of the form In^-^'^'^ j3. 



, and were first correctly obtained in [|T0|]. 
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3. Scale choices and theoretical uncertainties 

The starting point of the evolution of the soft function W/^, represented by the soft scale Hs 



appearing in (2.2), has to be chosen such that logarithms in the expansion in of sf are small, 
giving a stable perturbative behaviour at the low scale Us. For soft interactions, the natural scale is 
set by the kinetic energy of the tf pair, m^jS^. Accordingly, in the soft scale in the resummed 
cross section was set to 

Us = max[ksm,l5^^^,ksm,l5^] , (3. 1) 

with the constant chosen by default as kg = 2. In the upper interval, j8 > jScut, the soft logarithms 
InjS in the partonic cross section are correctly resummed by the running soft scale kstfifji^. If jScut 
is not too big, in the lower interval, j3 < j3cut, the frozen soft scale ksmt|5^^^^ still correctly resums 
the dominant contributions (~ InjScut) to the hadronic cross section, at the same time avoiding 
ambiguities related to the Landau pole in a^. A prescription for the choice of jScut is detailed in |^]. 

Contrary to the soft function, the hard function H^^, naturally lives at the parametrically larger 
scale given by the invariant mass of the tt pair. The default value for the hard scale /i/, appearing 
in the resummed hard function is therefore chosen as /i/, = 2mt. On the other hand, the scale 
of Coulomb interactions is set by the typical virtuality of potential gluons, ~ mj[5^. For the 



Coulomb scale jlc in (2.3) we thus choose 

/Xc = max[Cf asm?,2mij8] , (3.2) 

where the frozen scale at low values of j8, equal to the inverse Bohr radius CpCCsmt, signals the 
onset of tt bound-state effects. 

There is clearly some degree of ambiguity in the choice of the scales appearing in the re- 
summed result. Besides these, NNLL and approximated NNLO predictions are affected by uncer- 
tainties related to constant and power suppressed terms which are not controlled by resummation. 
Therefore, to reliably ascertain the residual theoretical error of the predictions presented below, we 
consider the following sources of ambiguity: 

Scale uncertainty: we vary all scales ^u, in the interval [/i,/2,2/i(] around their central values fii. 
He is varied while keeping the other scales fixed, /i/, and /X/ are allowed to vary simultane- 
ously, imposing the additional constraint 1 < < 4. For the fixed-order results (NLO 
and NNLOapp) the factorization scale /X/ and renormalization scale jj.^ are varied simultane- 
ously, with the constraint 1/2 < IXr/fXf < 2. The errors from varying {/Xy,/^/,} and jXc are 
added in quadrature. 

Resummation ambiguities: we consider three different sources of ambiguities: i) the difference 
between the default setting E = m,^^ compared to E = — 2mt in (pj|), ii) the difference 
between the NNLL implementation for the soft scale choices kg = \ ,4 to the default choice 
ks = 2, iii) the effect of varying /3cut by 20% around the default value for kg = 2 (see [||] for 
details). The resulting errors are added in quadrature. 



NNLO-constant: by default, the ^{cc^) constant in (^3]), C^^l ^, is set to zero. We estimate the 



pp , 

effect of a non- vanishing constant by considering variations Cp2 „ = ±|c|,2 with C^^l „ 

PP 5" PP )" PP 5^ 

the constants in the threshold expansion of the NLO cross sections. 
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Tevatron 


LHC =7 TeV) 


LHC (V^ =14 TeV) 


NLO 


^ ^0+0.36+0.51 


iro 1+18.5+13.9 
iJO.i_2i.2-13.1 


OO/1+107+65 
^^4_io6-58 


NNLOapp 


7 r.f^+0.21+0.69 
'•""-0.34-0.53 


161 1+12-3+15.2 


on 1+76+64 
^y^-69-63 


NNLL2 


7 99+0.31+0.71 
'•^^-0.47-0.55 


162 6+^-'^+^5-'^ 

^"^•"-7.5-14.7 


OQf,+40+65 
°^"-37-64 



Table 1: tt cross section at Tevatron and LHC from NLO, NNLOapp and NNLL2 approximations, for nit = 
173.3 GeV. The first eiTor set denotes the total theoretical uncertainty, the second the PDF+Us eiTor. 



In Section ^ the total theory error is obtained summing in quadrature the three aforementioned 
uncertainties. Additionally, we estimate the error due to uncertainties in the PDFs and the strong 
coupling, using the 90% confidence level set of the MSTW08NNLO PDFs and the five sets for 
variations of provided in [17]. This will be quoted separately from the theoretical error. 



4. Results 



We define our default prediction, NNLL2, by matching the NNLL resummed result to the 
approximated NNLO cross section, i.e. 

NNLL2 = NNLL - NNLL(a,'^) + NNLOapp , (4. 1) 

where NNLL(a^) is the expansion, up to order a^, of the resummed result. The top-quark mass is 
set to nit = 173.3GeV, and the renormalization and factorization scale are chosen as /ij = jj.,. = nif. 
The other scales are treated as explained in Section ^ For the convolution of the partonic cross 
section with the parton luminosity functions, Eq. ( |1.1[ ), we use the MSTW08 PDF set [17]. 

Numerical results for the cross section at the Tevatron and LHC for NLO, NNLOapp and the 
resummed NNLL2 implementation are given in Table [l], with the first error corresponding to the 
total theoretical uncertainty, computed as described in Section ^, and the second error the PDF+a^ 
uncertainty. The genuine NNLL corrections are sizeable and positive both at Tevatron (^^ +13%) 
and LHC (~ +9%), though their effect is partially compensated by switching from NLO to NNLO 
PDF sets for the NNLOapp and NNLL2 results, which results in a negative shift of the cross section. 
The bulk of the corrections beyond NLO is accounted for by the ^{(xj) terms, as evident from a 
comparison of NNLOapp and NNLL2. One can also notice a significant reduction of the theoretical 
uncertainty from NLO to NNL0app/NNLL2, both at Tevatron and the LHC, to the extent that for 
NNLOapp and NNLL2 the total error is dominated by the PDE+a^ error. 

The theory uncertainty of different approximations is shown in Figure [l], where the total 
theoretical error bands for NLO (dashed black), NNLOapp (dot-dashed blue) and NNLL2 (solid 
red) are plotted as functions of m,. At the LHC, one observes a nice convergence of the series 
NLO— >-NNLOapp— 5'NNLL2, with the resummed result having the smallest theoretical uncertainty, 
corresponding to about ±4.7%. At the Tevatron, however, the NNLOapp result exhibits the smallest 
error. We interpret this as an indication that the NNLOapp error is accidentally underestimated by 
the scale variation procedure. For NNLL2, the residual theoretical uncertainty is +4.3 — 6.5%. 

Predictions for the total cross section at the NNLL or approximated NNLO level have been 
published recently by several groups. These results differ in the resummation formalism adopted 
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Figure 1: NLO (dashed black), NNLOapp (dot-dashed blue) and NNLL2 (solid red) as function of ;«,. The 
bands correspond to the total theoretical uncertainty (that is, excluding the PDF+aj error) of the prediction. 

and in whether the total cross section is resummed, or differential distributions in different kine- 
matics limits, which are then integrated to obtain predictions for the inclusive cross section. They 
also differ in the treatment of Coulomb effects beyond NLO and of the constant terms at ff{aj), 
and some include sets of power-suppressed contributions in j8 . A comparison of the various results 
can thus give an estimate of the ambiguities inherent to resummation. This is shown in Figure 
where the predictions for NNLO and NNLL given here (black circles, [S]) are compared to the 



numbers by Kidonakis (blue up-pointing triangles, []130), Ahrens et al. (red diamonds, [ ]12p for 
NNLO and for NNLL in both IPI and PIM kinematics) and Cacciari et al. (green squares, 
[§]). The experimental measurements (purple down-pointing triangles, |[Tj, ^ ^, ^) are also given 
for comparison. At the LHC the predictions of different groups show a good agreement, with a 
somewhat better agreement for the NNLO results. At NNLL, the total spread of the four results 
is still smaller than the theoretical uncertainty of the NLO result. At Tevatron there appears to be 
a stronger tension between different results (both for central values and error estimates), with the 
envelope of all predictions at NNLO and NNLL being almost as large as the uncertainty of the 
NLO result. In ^ it was argued that this might be a consequence of the dominance, at Tevatron, 
of the qq production channel, which is less well approximated by its threshold expansion than the 
gg channel dominant at the LHC, thus leading to larger ambiguities in resummation. 

5. Mass determination 

Precise theoretical predictions of a,,- can be translated into measures of the top-quark pole mass 
nit. In first approximation, this can be done by comparing the mass dependence of the theoretical 
and experimental cross sections, and extracting the top mass at the intersection point. In a more 
sophisticated approach, one defines a likelihood function 

/(m,)oc J da fth{a\mt) ■ fe^p{a\mt) , (5.1) 

where /th and /exp represent normalized gaussian distributions centred around the theoretical pre- 
diction and measured experimental cross section, respectively, and extracts m, from the maximum 
of the probability distribution. A suitable parameterization of the experimental cross section in 
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cr[pbj 

9 r 



BFKS(»i, = 173.3 GcV) 
Kidoiiakis (m, = 173 GeV) 
Ahrens et al. (m, = 173. 1 GeV) 
Cacciari et al. (m, = 173.3 GeV) 



ft 

t 



<r[pbl 
220 r . 



BFKS(»i, = 173.3 GeV) 
Kidonakis (m, = 173 GcV) 
- ♦ Ahrens et al. (m,= 173.1 GeV) 
Cacciari et al. (m, = 173.3 GcV) 



LHC (7 TeV) 



I}* *1 



Figure 2: Comparison of different NNLO and NNLL predictions, see the text for explanation and refer- 
ences. The error bands include theoretical uncertainties, but no PDF+aj errors. The rightmost set of points 
represents the most recent experimental measurements, which assume a top mass of 172.5 GeV. 



o-„(m,)[pb] 




Figure 3: mt dependence of the experimental tt cross section [||] (solid black) and of the NNLL2 approxi- 
mation presented here (solid red). The dashed lines represent the total uncertainties of the two curves. 

terms of has been provided by the ATLAS collaboration and is plotted in Figure ^. Using 
the NNLL2 prediction presented here as theoretical input in (5J^) and the data from gives the 
following value for the top pole mass, 

mt = m.ittjGey, (5.2) 
which agrees with the direct-reconstruction measurement of Tevatron, nit = 173.3 ±1.1 GeV. Note 



that the error in ( |5.2[ ) does not include the uncertainty deriving from identifying the Monte Carlo 
mass in the experimental input with the pole mass nit. Allowing for a ±2 GeV difference would 
correspond to an additional uncertainty of ±0.65 GeV on the extracted mass. 

6. Conclusions 

We presented new predictions for the tf total production cross section which include NNLL 
resummation of soft and Coulomb effects. Resummation increases the cross section relative to the 
fixed-order NLO result, and leads to a significant reduction of the theoretical uncertainty. Our 
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numbers are in good agreement with experimental measurements at the Tevatron and LHC. At the 
LHC they also show a reasonably good agreement with other NNLL predictions, though bigger 
differences are found at the Tevatron. 

References 

[1] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 105 (2010) 012001 [arXiv: 1004.3224 
[hep-ex]]; CDF Collaboration, conference note 9913, 2009. 

[2] V. M. Abazov etal. [DO Collaboration], Phys. Rev. D84 (2011) 012008 [arXiv: 1101. 0124 [hep-ex]]; 
V. M. Abazov etal. [DO Collaboration], Phys. Lett. B704 (2011) 403 [arXiv:1105.5384 [hep-ex]]. 

[3] G. Aad et al. [ATLAS Collaboration], arXiv: 1108.3699 [hep-ex]; ATLAS Collaboration, conference 
note ATLAS -CONF-20 11 - 1 2 1 , 20 11 . 

[4] S. Chatrchyan et al. [CMS Collaboration], Phys. Rev D84 (201 1) 092004 [arXiv: 1 108.3773 
[hep-ex]]; CMS Collaboration, conference note CMS PAS TOP-1 1-024 

[5] E. Laenen, J. Smith, and W. L. van Neerven, Nucl. Phys. B369 (1992) 543-599; S. Catani, M. L. 
Mangano, P Nason, andL. Trentadue, Phys. Lett B378 (1996) 329-336, [arXiv:hep-ph/9602208]; 
E. L. Berger and H. Contopanagos, Phys. Rev D54 (1996) 3085-31 13, [arXiv:hep-ph/9603326]; 
N. Kidonakis, J. Smith, and R. Vogt, Phys. Rev D56 (1997) 1553-1570, [arXiv:hep-ph/9608343]; 
R. Bonciani, S. Catani, M. L. Mangano, and P Nason, Nucl. Phys. B529 (1998) 424^50, 
farXiv:hep-ph/9801375]; N. Kidonakis, E. Laenen, S. Moch, and R. Vogt, Phys. Rev D64 (2001) 
114001, [arXiv:hep-ph/0105041]. 

[6] T. Becher and M. Neubert, Phys. Rev D79 (2009) 125004 [Erratum-ibid. D80 (2009) 109901] 
[arXiv:0904.1021 [hep-ph]]; M. Beneke, P Falgari and C. Schwinn, Nucl. Phys. B828 (2010) 69 
[arXiv:0907.1443 [hep-ph]]; M. Czakon, A. Mitov and G. F Sterman, Phys. Rev D80 (2009) 074017 
[arXiv:0907. 1790 [hep-ph]]; A. Ferroglia, M. Neubert, B. D. Pecjak and L. L. Yang, Phys. Rev. Lett. 

103 (2009)201601 [arXiv:0907.4791 [hep-ph]]. 

[7] V. Ahrens, A. Ferroglia, M. Neubert, B. D. Pecjak and L. L. Yang, arXiv: 1103.0550 [hep-ph]. 

[8] M. Beneke, R Falgari, S. Klein and C. Schwinn, Nucl. Phys. B855 (2012) 695 [arXiv: 1109. 1536 
[hep-ph]]. 

[9] M. Cacciari, M. Czakon, M. L. Mangano, A. Mitov and P. Nason, arXiv:l 11 1.5869 [hep-ph]. 

[10] M. Beneke, M. Czakon, P Falgari, A. Mitov and C. Schwinn, Phys. Lett. B690 (2010) 483 
[arXiv:0911.5166 [hep-ph]]. 

[11] M. Aliev, H. Lacker, U. Langenfeld, S. Moch, P. Uwer andM. Wiedermann, Comput. Phys. Commun. 
182 (2011) 1034 [arXiv:1007.1327 [hep-ph]]. 

[12] V. Ahrens, M. Neubert, B. D. Pecjak, A. Ferrogha and L. L. Yang, Phys. Lett. B703 (201 1) 135 
[arXiv:1105.5824 [hep-ph]]. 

[13] N. Kidonakis, arXiv: 1109.3231 [hep-ph]. 

[14] M. Beneke, P Falgari and C. Schwinn, Nucl. Phys. B842 (2011) 414 [arXiv: 1007.5414 [hep-ph]]. 

[15] T. Becher, M. Neubert and G. Xu, JHEP 0807 (2008) 030 [arXiv:0710.0680 [hep-ph]]. 

[16] M. Czakon and A. Mitov, Nucl. Phys. B 824 (2010) 111 [arXiv:0811.4119 [hep-ph]]. 

[17] A. D. Martin, W. J. Stirling, R. S. Thorne, and G. Watt, Eur. Phys. J. C63 (2009) 189-285, 

[arXiv:0901.0002 [hep-ph]]; A. D. Martin, W. J. Stirling, R. S. Thorne, and G. Watt, Eur Phys. J. C64 
(2009) 653-680, [arXiv:0905.3531 [hep-ph]]. 



8 



